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Devising a statistical mechanics framework for jammed granular materials is a challenging task as 
those systems do not share some important properties required to characterize them with statistical 
thermodynamics tools. In a recent paper [Asenjo et al. PRL, 2014], a new definition of a granular 
entropy, that puts the protocol used to generate the packings at its roots, has been proposed. Fol¬ 
lowing up these results, it is shown that the protocol used in [Asenjo et al. PRL, 2014] can be recast 
as a canonical ensemble with a particular value of the temperature. Signature of gaussianity for 
large system sizes strongly suggests an asymptotic equivalence with a corresponding microcanonical 
ensemble where jammed states with certain basin volumes are sampled uniformly. We argue that 
this microcanonical ensemble is not Edwards’ microcanonical ensemble and generalize this argument 
to other protocols. 

PACS numbers: 45.70.-n,45.70.Cc,05.90.-|-m 


I. INTRODUCTION 

Granular materials are large assemblies of particles 
which are athermal and dissipative and as such their sta¬ 
tistical properties can hardly be characterized by using 
the standard tool box of statistical thermodynamics. It 
has then been considered a challenge in the community 
to devise a novel statistical framework for these systems 
that would have both explanatory and predictive powers 
W- One of the most appealing routes towards a sta¬ 
tistical framework to characterize jammed states is the 
one proposed by Sam Edwards more than two decades 
ago which hypothesizes that if one has a protocol to gen¬ 
erate (granular) packings of N particles at some fixed 
volume V or packing fraction (j), then each jammed pack¬ 
ing will be as likely as any other to appear with proba¬ 
bility ~ where fljammed{N,(j)) is the to¬ 

tal number of mechanically stable structures that can 
be made out of these N particles confined in the pre¬ 
scribed volume V — EE]. Although such a claim does 
not bear any micro-mechanical justification yet, there is 
still a possibility that, akin to throwing a six facet dice 
(system for which we do not have any micro-mechanical 
justification either), it so happens that most protocols 
could well be described by a probability measure which 
assumes a uniform sampling of the states in virtue e.g. 
of Pascal’s principle of indifference. Still, Edwards’ pro¬ 
posal contrasts deeply with studies like that of Jiao et al. 
[315] which showed that one could tailor many packing¬ 
generating protocols that lead to any possible packing 
fraction corresponding to various degrees of structural 
randomness for the obtained jammed states (hence rul¬ 
ing out e.g. packings with the same packing fraction but 
with a different randomness). Also, recent advances on 
the behaviour of hard sphere glasses have permitted to 
probe the glass phase close to the jamming transition. 
This probing can be done in two different ways: either 
by focusing on the “ typical ” (in a statistical thermody¬ 
namic sense) states that can exist close to jamming |9] 
or by following initial states from a dilute phase m up 


to jamming m- Interestingly, these two different ways 
of choosing glassy states yield a different physics which 
strengthen the idea of what seems to be a protocol depen¬ 
dent measure. It is therefore important to devise tests to 
assess whether Edwards’ conjecture is rather the rule or 
the exception to the rule. To this aim, different strate¬ 
gies have been used. On the one hand, instead of testing 
Edwards’ microcanonical hypothesis, it has been more 
practical to test a canonical extension of it which ought 
to apply to packing-generating protocols that do not fix 
exactly the volume occupied by the jammed states. In 
this case however the various theoretical [IM3 , numer¬ 
ical UMI] and experimental [501 155H55] studies done on 
the subject were not, overall, very conclusive. That is 
because different protocols and assumptions were used 
to apply Edwards’ canonical ensemble to practical cases. 
On another hand, it has been proposed for a long time 
to use soft materials as models for granular materials 
[5H150055] . In this case, the idea consists in interpreting 
the discrete set of jammed granular states as being min¬ 
ima of a potential energy surface (also referred to as an 
energy landscape; term that we should use interchange¬ 
ably with potential energy surface throughout the paper) 
of soft particles [33]. This allows then both experimental 
[501155] and numerical [55] [50] [55] [50] realizations of soft 
material models of granular packings that can be used in 
principle to test Edwards’ original hypothesis. 

Now, it appears that even the first set of studies look¬ 
ing at granular jammed states with a non-constrained 
packing fraction can be interpreted in terms of a land¬ 
scape that would be a volume or density landscape m 
where the minima of this landscape would be the set of 
states for which, under a prescribed dynamical model 
that prevents particle overlaps {e.g. the Lubachevsky- 
Stillinger (LS) algorithm [58] [39]), it is impossible to 
reach a lower volume configuration before first increasing 
the volume accessible to the particles [ID]. This allows 
then a rather unified picture of two of the most ubiq¬ 
uitous ways to generate jammed states expressed in the 
vocabulary of landscape exploration. Figd] gives an il- 
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Figure 1. (Color online) Schematic representation (dotted 
line) of an energy or volume landscape that can be used to 
generate four mechanically stable packings (A, B, C, D). For 
each minimum, it is not possible to decrease further the cost 
function without first increasing the volume available (volume 
landscape interpretation like the LS algorithm) or the energy 
of the system (energy landscape interpretation). The solid 
blue line represents the basin of attraction of the minimum C 
while the length of the red segment B'C' represents its size. In 
a multidimensional space, this size becomes a hyper-volume 
in the space of coordinates. 


lustration of the idea underlying this interpretation of 
jammed states as being either volume minima in a den¬ 
sity landscape of hard particles or energy minima of the 
potential energy surface of a soft material model of gran¬ 
ular matter. 


In this article, we will mostly focus on the energy land¬ 
scape (potential energy surface) exploration as we seek 
physical consequences — regarding e.g. Edwards’ con¬ 
jecture — of numerical results obtained in [3S]. To this 
end, we hrst recall very standard results on the equiva¬ 
lence between the canonical and the microcanonical en¬ 
sembles in the thermodynamic limit. We then move to 
jammed systems as studied in |35] . After having recalled 
the protocol they used and the results they obtained, we 
will show that there is a strong analogy between an in¬ 
finite number of energy landscape-based protocols and 
a — granular — canonical ensemble. We will moreover 
show that, akin to what happens in statistical thermo¬ 
dynamics, this canonical ensemble tends toward a micro- 
canonical measure. We discuss the implications of this 
equivalence and claim that this microcanonical measure, 
in general, does not coincide with the Edwards’ measure. 


II. EQUIVALENCE BETWEEN THE 
MICROCANONICAL AND CANONICAL 
ENSEMBLES 

This section is dedicated to recall how does the canon¬ 
ical ensemble become equivalent to a microcanonical 
ensemble in the thermodynamic limit. Although, the 
derivation that follows can be found in most text books 
on statistical thermodynamics {e.g. in mi), we briefly re¬ 
derive it here to introduce notations and concepts that 
will be useful in the remaining of the manuscript. We 
start off by writing the partition function in the canoni¬ 
cal ensemble, for N particles in a volume V with inverse 
temperature j3 = l/ksT, via a sum over energies E i.e.: 

1 r+oo 

Q{N, = dE cu{E, V, N)e-^^ (1) 

where uj{E, V, N)dE is the total number of states in 
the energy range [E,E + dE]. The idea consists then 
in defining the microcanonical entropy S{E, N,V) via 
u}{E) = and it comes that the partition 

function reads: 


^+oo 

Q{N,V,P)= / dE (2) 

Jo 

where A{E) = E — TS{E) + fcsTln V! is the Helmholtz 
free energy of the system at energy E and where the 
-l-lnV! term ensures extensivity of A. To formally ad¬ 
dress the extensivity property, we introduce the notion 
that two functions f{N,V,f3) and g{N,V,f3) are equiv¬ 
alent in the thermodynamic limit (and noted f g) 
iff limjv,v->.-i-oo g ~ ^ fixed packing fraction cj) = 

Nvq/V ~ 0{\) being a length scale related to 

the size of the particles; where d stands for the spa¬ 
tial dimension). A being extensive, then — by defini¬ 
tion — there exists a function a{e,(p) ~ 0(1) such that 
A{E) ~ Na{e, (j)) and where e and a are the energy and 
free energy per particle. In the thermodynamic limit, Eq. 
([^ can then be recast as: 

r+oo 

Q{N,V,p)^N (3) 

Jo 


Now, if a has a minimum at say e*, then this integral 
can be evaluated by using a saddle point approximation. 
Assuming this extremum lies on the real line, the idea is 
to expand a(e) = a(e*) -|- ^a"{e*){e — e*)‘^ + 0{{e — e*)^). 
Replacing it in Eq. up to the second order gives: 

Q(V,U,^) ~ r de 

Jo 

From Eq. Q, we see that all the statistics is now es¬ 
sentially in a gaussian form. From very standard statis¬ 
tical mechanics, we know that Var(E) oc Cy where Cy 
is the specific heat capacity at fixed volume. Now, the 
heat capacity is also an extensive quantity and therefore 
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Var{E) ~ N. Now, in the case of Eq. Q, the statis¬ 
tics is that of e and not E. This is no issue as we can 
use the fact that Var{e) = Var{E)/N'^ ^ 1/N. It then 
implies, in case we had any doubts, that a"(e*) ~ 0(1)- 
Evaluating the gaussian integral then yields: 


(5) 

where K ~ 0{\) is a constant. Hence, in the thermo¬ 
dynamic limit, one finds that A{N,V,f3) ~ Na{e*) = 
N{e* — Ts{e*)) (where s contains here part of the fV! 
introduced in the definition of A in Eq. <§)■ Everything 
is thus as if the energy was fixed at the value Ne*, from 
which we subtract the corresponding (extensive) micro- 
canonical entropy. The role of the thermal bath is en¬ 
capsulated in the possibly complicated relationship be¬ 
tween e* and /?. Now, since the partition function is the 
generating functional of all the moments of the canon¬ 
ical distribution, we have therefore retrieved that there 
is a statistical equivalence between the canonical ensem¬ 
ble and a microcanonical one in the thermodynamic limit 
i. e. there exists a one to one correspondence between the 
physics going on at some /3 value in the canonical ensem¬ 
ble and the physics characterized by the corresponding 
microcanonical ensemble aX E = Ne*{f}). 

Another way to see this equivalence is to look at the 
integrand in Eq. and see that we can build a prob¬ 
ability distribution for e from it by simply dividing the 
gaussian weight by the partition function itself. Since we 
have calculated it in Eq. ([^, we see that at least up to 
order 2, we have that: 


PAr(e) 



KN{e-e*f\ 

2 ) 


( 6 ) 


One can easily check that pat is essentially a sequence of 
functions whose limit is the Dirac delta distribution such 
that: 


lim pAr(e) = 6{e — e*) (7) 

which is essentially a microcanonical measure at hxed 
energy Ne*. 

III. GENERATING JAMMED STATES AND 
THE LANDSCAPE AMBIGUITY 

The problem of generating model granular jammed 
states from an energy landscape can be divided into 
two parts. The first one consists in identifying a model 
whose energy minima mostly correspond to — granular 
— jammed states. The second one consists in providing a 
procedure for the landscape exploration that ultimately 
decides the frequency with which each jammed state is 




(a) 








Figure 2. Jammed state as generated by potential energy 
minimization. One start by picking a state (a) from an equi¬ 
librated hard-sphere fluid. Then, the soft shell interaction 
is turned-on and the packing fraction is now such that there 
are many overlaps (b). The system relaxes then to the cor¬ 
responding potential energy minimum, yielding the jammed 
state (c). 


visited. For soft matter models, the potential energy sur¬ 
face is that of a system of hard core particles — whose 
radii can go to zero in some models — covered by a soft 
shell — using hertzian [551 HIl SH I13]j elastic [23] or 
WCA [35] potentials —. For our purposes, we recquire 
that the packing fraction of the whole system, when de¬ 
fined with respect to the total particle sizes {i.e. hard 
core + soft shell), is way above jamming. The “ above 
jamming ” condition ensures that all minima correspond 
to jammed states and not to a part of the conflguration 
space that would be fluid-like. That way, the jammed 
basins tile certainly the underlying fluid space; which is 
a central feature of the method as we shall see below. 
The price to pay is then that the present analysis can¬ 
not straightforwardly be used to probe the whole energy 
landscape at the jamming transition. Study of this par¬ 
ticular point below which most disordered hard particle 
assemblies are fluid [H] goes beyond the scope of the 
present paper and is still the object of intensive studies 
P [211261 [311132 Sa ESSS] . 

Now, regarding the exploration procedure, it is very com¬ 
mon to choose a packing fraction for the hard core part 
of the system that is sufficiently low to avoid any glassy 
behaviour. In that case, if one picks, at random, config¬ 
uration states from an equilibrated hard core fluid, then 
this is equivalent to exploring uniformly the potential 
energy surface associated to the soft shell interactions. 
This protocol is illustrated with typical simulation steps 
in Fig. H 

A consequence of this exploration protocol is that, 
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upon a quench to the “ nearest ” minimum, a fluid state 
picked in this way is more likely to be captured by a 
minimum with a big basin of attraction rather than by 
a small one. In fact, the probability to end up in a par¬ 
ticular basin of hyper-volume V is exactly V/Vfluid with 
this aforementioned protocol; where Vfluid is the volume 
tiled by all the basins. Thus, in addition to being simple, 
this widely used protocol is also very convenient because 
we can infer straightforwardly what the jammed states 
statistics has to be (and it is not uniform). 

Now, to prevent any confusion, a comment is here in 
order. The energy landscape description of soft material 
models of a granular medium discussed above (and in 
FigQ looks very similar to the free energy landscape de¬ 
scription of structural glasses. In fact, many studies have 
looked deep into the glassy regime with the use of the 
Replica Symmetry Breaking tools to probe the jamming 
transition for soft and hard particles {e.g. [2l 1^ IdSUhO] !. 
In particular, it has been found recently that when us¬ 
ing the smallest resolution possible to define the basins 
of this free energy landscape for hard spheres, the actual 
hyper-volume of these basins was vanishing at jamming 
[5]. This contrasts with the claim of the previous para¬ 
graph where the hyper-volume of the basins of attraction 
of jammed states have no reason to be zero (this is quite 
clear in Figtoo). These two observations are in fact 
non-contradictory for they do not talk about the same 
landscapes. On the one hand, the free energy landscape 
of structural glasses — which is a Gibbs entropy land¬ 
scape for hard spheres — is an emerging property of the 
thermodynamic state of the system; that is, upon fol¬ 
lowing a particular branch of a metastable fluid phase, 
some regions of the configuration space become more and 
more dynamically disconnected from one another as we 
go deeper into the glass phase. This gives rise to the suc¬ 
cessive appearance of more and more “ structural basins 
” which characterize these effectively disconnected re¬ 
gions within which, in the glassy regime, particles are 
still free to move and hence a thermodynamic Gibbs en¬ 
tropy can be associated to them. It is then no surprise 
that, when approaching jamming from the glass, the “ 
typical ” entropy per basin vanishes because by defini¬ 
tion, no particle has any freedom to jiggle at jamming 
and the Gibbs partition function has to be zero there. 
That is simply because, at jamming, the dimensional¬ 
ity of the space to explore decreases — essentially from 
continuous to discrete [33] — and as a consequence, the 
Gibbs measure in the continuous configuration space of 
hard spheres has to be exactly zero for any set in this new 
space of jammed states. On another hand, the potential 
energy surface discussed in Ref. [3S] and in the present 
paper is a feature of a soft material model of a granular 
material and hence of the effective protocol used to gen¬ 
erate model granular jammed states. The definition of 
basins of attraction — and their hyper-volume — then 
differs greatly in these two cases. The free energy land¬ 
scape of glasses would define a basin as being the set of 
all fluid configurations compatible with a particular aver¬ 


age structure (set that ought to be of cardinality one at 
jamming) whose hyper-volume is then often related to a 
typical cage size reminiscent of cell model methods used 
to compute the Helmholtz free energy of crystals. Essen¬ 
tially, the smaller the cage, the smaller the basin entropy 
and the harder it becomes to compress the system any 
further. Upon following a metastable branch of the fluid 
phase, this leads to a divergence of the thermodynamic 
pressure at a packing fraction much lower than that of 
close packing which is symptomatic of a vanishing typ¬ 
ical cage size as discussed by Kamien and Liu |33|. A 
contrario, the hyper-volume of a basin of attraction of 
a potential energy surface would be the set of all initial 
states which lead, upon a prescribed quenching dynam¬ 
ics, to the same jammed structure and has absolutely no 
reason to be of zero measure. For a glass, it could be the 
set of initial fluid states at an initial equilibrium pack¬ 
ing fraction (fi which, upon quenching to a higher final 
packing fraction (j)f beyond the glass transition, reach the 
same average structure. This is very different from the 
number of fluid states at (j>f which are compatible with 
some average structure. Now, in our case, looking at Fig. 

the basin of attraction of the jammed structure (c) is 
the set of all (a)-states leading to the structure (c) upon 
following the steps decribed in the above paragraph and 
summarized in Fig. This discussion on the difference 
between basin entropies in a glassy free energy landscape 
close to jamming and basin hyper-volumes of a potential 
energy surface echoes a more general discussion on other 
— non-thermodynamic — entropies that can be defined 
to characterize jammed states m and had to be made 
clear before moving on to the main point of the paper. 


IV. CANONICAL ENSEMBLE FOR JAMMED 
STATES 


Ref. [35] aimed at estimating the total number of 
jammed states for a given packing fraction (f and num¬ 
ber of particles N (for a poly-disperse system). To do 
so, the authors imagine an ideal — unknown — protocol 
which would sample all of the jammed states uniformly 
and define the mean basin hyper-volume as being: 


(V)^ 




led 


INA) 




^jammed 


( 8 ) 


The idea is then to infer an unbiased distribution Pu(,F) 
for the free energy F = — ln(VA^!) [52] of a basin such 
that: 


(V) = J dFpu{F)e-^ (9) 

The major problem is that, as specified before, Pu{F) is 
unknown since the protocol described above to actually 
generate the jammed states is far from sampling them 
uniformly. At this stage, it is worth pointing out that 
the ideal unbiased algorithm that would sample all the 
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jammed states uniformly has a priori nothing to do with 
the basin size distribution p„(F). Indeed, using an energy 
landscape is just a particular way — among others — 
to get jammed states. So, while we can always imagine 
jammed states as being minima of an energy landscape, 
we could have used another feature of jammed states to 
sample them uniformly, and that would not have changed 
their basin volume distribution. Let us now denote Pb{F) 
the free energy distribution of the basins as sampled by 
the protocol with a bias proportional to the hyper-volume 
of the basins. It is clear that: 


Pb{F) cx Pu{.F)V{F) = Pu{F)e-^, (10) 


and now something interesting happens. The expres¬ 
sion in Eq. 0 can indeed be interpreted as a partition 
function similar to that of Eq. 0 associated to the bi¬ 
ased weight in Eq.(lO) with which each basin is sampled 
with the protocol in use. The interpretation goes then 
as follows: everything is as if we had “ microstates ” 
which are the jammed states and for each of them we 
can measure an “ energy ” which is nothing but the free 
energy associated to their corresponding basin of attrac¬ 
tion {F = — ln(VA^!)). In Eq. p„(F) plays exactly 
the same role as uj{E) in Eq. ([^ corrected by a Boltz¬ 
mann weight penalizing high energies in the case of Eq. 
Q and small volumes in the case of Eq. (10). In addi¬ 
tion, F is here strictly bounded from below because of the 
size of the box but has no obvious strict upper bound. A 
major difference with usual statistical mechanics though, 
is that both the biased and the unbiased distributions for 
F have to be integrable and in general, we assume that 
most of their moments — if not all — have to be defined. 
A consequence of these general requirements is that the 
free energy density of state noted ujjammediF) oc p„(F) 
has to be a fast decaying function for both small and 
large values of F (almost defining a compact support). 


Finally, for the analogy with a canonical ensemble to 
be complete, we can wonder what is the equivalent of 
/? in the case of jammed states (that we denote (do)- It 
turns out that Pc is related to how sensitive the sampling 
protocol is to the size of the basins. With the protocol 
described above, the sampling selects proportionally to 
the hyper-volume of the basin and /3 g = 1, but one could 
imagine other types of sampling for which the bias with 
respect to uniformity would be some power of the hyper¬ 
volume of those basins and this exponent would be the 
corresponding Pq. 


V. TOWARDS A MICROCANONICAL 
ENSEMBLE 

One of the results of Ref. [3S] is that the biased distri¬ 
bution Pb{F) tends towards a gaussian distribution for 
large system sizes. It is also found that {F)b ~ Nf* i.e., 
the biased average of F increases linearly with the sys¬ 
tem size [3S]. In addition, it is found that Var{F) ~ Nc 


(where c ^ 0{1)) so that in the end, if we look at the bi¬ 
ased distribution of the free energy per particle / = F/N, 
it comes (for large system sizes): 


Pbif) 






( 11 ) 


As we have seen before, such a function will tend to¬ 
wards a Dirac delta distribution and, in effect, will be 
uniformly sampling basins with a single free energy that 
is Nf*. Now, the corresponding ensemble is one where 
the free energy is essentially Nf* for every basin sam¬ 
pled. It is worth stressing that such a claim was al¬ 
ready hypothesised by Wang et al. [23] to justify that 
Edwards’ measure would be retrieved in the thermody¬ 
namic limit even for biased sampling protocols. Although 
we acknowledge that the shrinking of some basin-related 
probability distribution was indeed inspired, we disagree 
with the conclusion that the Edwards measure will al¬ 
ways be retrieved in the thermodynamic limit and shall 
explain why below. 

There are in fact two ways to interpret Eq. (1^: a) 
all possible jammed states are sampled uniformly because 
the landscape becomes such that all basins have the same 
hyper-volume corresponding to the free energy N f* or b) 
only a subset of all possible states is sampled uniformly in 
the thermodynamic limit and that is the one comprising 
basins with a free energy Nf*. The remaining part of 
the paper will deal with this issue. 


VI. FINDING THE DENSITY OF STATE UPON 
UNBIASING Pi 


Let us first remark that, if the biased distribution is 
genuinely gaussian (which is what we will assume here 
based on Ref. [35]), then so is the unbiased one. To 
see this, one simply needs to use Eq.(lO) together with 
the following identity —a{x — x)^ + x = —a(x — (x + 
(2a)~^))^ — X — (4a)“^. As a matter of fact, since we 
multiply Pb(,F) by an exponential to retrieve p„(F), its 
mean will be shifted towards higher values of the free en¬ 
ergy (an illustration of this is given e.g. in Ref. m)- 
Now, the fact that the density of state of F (or equiv¬ 
alently the probability distribution pu{F)) tends to be 
gaussian for large system sizes was also found in [33| by 
sampling all the minima of small binary mixtures of elas¬ 
tic disks. This suggests that the gaussianity of Pu{F) 
(or UjjammediF)) does not depend very much on the ac¬ 
tual energy landscape as long as the interactions are kept 
short ranged and isotropic. 

In the end, for such energy landscapes, the free en¬ 
ergy per particle distribution — as measured by uniform 
sampling — is a delta function centered around f*+A 
where A accounts for this shift and is a positive number 
unimportant for the present discussion. 

We already see here that if we believe in the interpre¬ 
tation a) above, that is, interpreting Eq. 0 as meaning 
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that all states are sampled uniformly and have free en¬ 
ergy Nf*, then there is a contradiction with the fact that, 
when all the jammed states are sampled uniformly and 
we infer the corresponding unbiased distribution from 
Pbi then we know that there are overwhelmingly more 
states with free energy N{f* -\- A) than with free energy 
Nf*. Thus, interpretation a) cannot be right. The only 
remaining choice, common from a statistical mechanics 
point of view, is that to a canonical sampling with /?g = 1 
corresponds an equivalent uniform sampling of a partic¬ 
ular subset of all the jammed states that contains only 
jammed states with free energy Nf*. 

As said before, one can imagine protocols whose sen¬ 
sitivities to the basin size will go as . In such a case, 
the problem is equivalent to that of a canonical ensem¬ 
ble with an inverse “ temperature ” Pc and a general 
partition function: 

Qg{N,4>,Pg) = j dF iO,ammed{F)e-^^^ ( 12 ) 


The crucial point to make here is that the gaussianity of 
t^jammed{F) (or equivalently Pu{F)) for large N is not 
protocol-dependent and is a feature of the soft matter 
model. As a consequence, any canonical weight with a 
temperature /3g, in virtue of the trivial extension of Eq. 
(10) to any Pg value, will tend towards a gaussian for 
large system sizes and will be statistically equivalent to a 
uniform measure within a corresponding free energy slice 
in the thermodynamic limit. 


In contrast to what was suggested in Ref. [23], we 
stress that, in general, the limiting microcanonical distri¬ 
bution discussed above does not coincide with Edwards’ 
ensemble as the latter requires that all possible states at a 
given N and (j) are uniformly sampled. In fact, Edwards’ 
measure would correspond to the sum of u}{F) over all 
basin free energy slices i.e. over all the possible limiting 
microcanonical ensembles as derived in the present pa¬ 
per. Now, an interesting case is that of Pg = 0. This 
corresponds to an infinite temperature in our interpreta¬ 
tion of the protocol and indeed in this case there is no 
bias and all states can be sampled uniformly; the Ed¬ 
wards measure is then retrieved as a very specific case of 
our interpretation. That being said, it is worth mention¬ 
ing that an interesting discussion on how fast a quench 
is performed in an energy landscape and its effect on the 
statistical mechanics description of granular media has 
been proposed in [24|. It was indeed argued that slowly 
quenched protocols were compatible with Edwards’ mea¬ 
sure while, at least for small systems, fast quenches were 
not. This claim seems quite un-intuitive at first sight 
since we know that e.g. for hard spheres, a slow quench in 
packing fraction simply leads to an entropy-driven crys¬ 
tallization that ends up in a jammed close-packed ECC 
or HCP structure; this is not quite what one would call 
a uniform sampling of jammed states. This intuition 
is indeed confirmed by running standard Monte Carlo 
and Molecular Dynamics simulations but also LS sim¬ 


ulations which, at sufficiently small compression rates, 
can mimic equilibrium behaviour laiss]. Now, the nu¬ 
merical quenching protocol of Ref. [231 also of Refs. 
[271 [53] with similar findings) actually differs from a usual 
LS algorithm in that it allows each particle to be sub¬ 
ject to a background fluid-like dissipation that forces the 
system to eventually reach a mechanical equilibrium — 
not necessarily jammed — after each compression step 
[54] . Motion of the particles is then only triggered by 
particle overlaps upon further compression. Such a pro¬ 
tocol is very similar in spirit to the one of Ref. [3T] cau¬ 
tiously designed to study the jamming of soft athermal 
systems. Hence, it is slightly misleading to interpret such 
a protocol as being related to a single potential energy 
surface — that of the final packing fraction — for the 
system explores many different energy landscapes along 
its compression trajectory. This phenomenology of ever- 
changing landscape upon reaching jamming bears some 
resemblance with what happens when trying to probe 
jamming from model spin glasses m or more recently 
hard sphere glasses [9], although, because the particles 
remain forever soft in the studies |24l I54j . a non-zero 
hyper-volume of the basins of attraction is expected (and 
indeed measured). 

In the present formalism, we claim that slowly reach¬ 
ing a minimum of an energy landscape (in the sense dis¬ 
cussed in the previous paragraph) from an initial fluid 
configuration may in fact well be an artificial way of de¬ 
creasing the value of (3g towards zero but it is however a 
bit too hasty to conclude that the corresponding statis¬ 
tical measure is that of Edwards’. Instead, we propose 
that, as the protocol used is more and more able to ef¬ 
fectively decrease /3 g, one should observe more and more 
compatibility with Edwards’ proposal and reach equality 
only in the limiting case /3 g = 0. 

Einally, in case one were to be surprised that we can 
find an equivalent microcanonical measure — centered 
at N{f* -I- A) — even when the sampling is unbiased 
and probing basins of all sizes i.e. when /?g = 0, it 
is worth recalling that the same goes for any thermody¬ 
namical system with a bounded energy spectrum in usual 
statistical thermodynamics. As a prototypical example, 
in a paramagnetic spin system coupled to an external 
magnetic field at inverse thermodynamic temperature /I, 
one easily finds that the canonical ensemble at /3 = 0 
is equivalent to a microcaninical ensemble where the en¬ 
ergy E = (E)^=o = 0. In this spin example, and in our 
granular case too, the unbiased sampling over all pos¬ 
sible states is overwhelmingly dominated by the states 
with the most probable energy (or by basins with the 
most probable free energy in our granular case). 


VII. CONCLUSION 

In this article, we have shown that it was possible to 
interpret the role of a widely used numerical protocol to 
generate jammed states as a particular “ granular” tern- 
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perature value — /3 g = 1 — of a much richer canonical 
ensemble. We have also shown that, assuming results 
from [5^ can be extended to more general cases of en¬ 
ergy landscape based protocols (cf. those used in Ref. 
[24j 1. this canonical ensemble is equivalent to a micro- 
canonical ensemble in the thermodynamic limit where 
some jammed states are then sampled uniformly. We ar¬ 
gued that Edwards’ microcanonical ensemble was a very 
specific case of this richer picture where the protocol in¬ 
verse temperature is /?g = Oj otherwise the other granu¬ 
lar microcanonical measures Pg > 0 do not coincide with 
Edwards’. 

These results shed light on how to take the role of the 
protocol seriously and yet be able to formulate the prob¬ 
lem in a formalism identical to that of statistical ther¬ 
modynamics. We hypothesise that the results discussed 
in this paper for energy landscape based protocols could 
extend to volume landscapes m and would give a new 
intuition about the statistical mechanics of tapping pro¬ 
tocols for instance. It is worth noting that, even in the 
case of a volume landscape, the probability assigned to 


each jammed state would be directly linked to the hyper¬ 
volume of the corresponding basin of attraction. Only if 
it is possible to have a mapping from these hyper-volumes 
of basins of attraction to the corresponding packing frac¬ 
tions of the jammed states, can we try to infer from it 
some statistical predictions about the probability distri¬ 
bution of volume or packing fraction for a given protocol. 
Finally, the interpretation here introduced opens the win¬ 
dow for novel granular simulations in which an energy 
landscape exploration could be coupled to a Metropo¬ 
lis criterion enforcing the value of Pq. This could give 
much more freedom in terms of protocols used to gener¬ 
ate jammed states and interpret their physical properties. 
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